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ABSTRACT 

We present an asteroseismic study of the solar-like stars KIC 11395018, KIC 10273246, KIC 10920273, KIC 10339342, and 
1 J - KIC 11234888 using short-cadence time series of more than eight months from the Kepler satellite. For four of these stars, we 

' derive atmospheric parameters from spectra acquired with the Nordic Optical Telescope. The global seismic quantities (average large 

JH , frequency separation and frequency of maximum power), combined with the atmospheric parameters, yield the mean density and 

surface gravity with precisions of 2% and ~0.03 dex, respectively. We also determine the radius, mass, and age with precisions of 
2-5%, 7-11%, and ~35%, respectively, using grid-based analyses. Coupling the stellar parameters with photometric data yields an 
asteroseismic distance with a precision better than 10%. A vsin; measurement provides a rotational period-inclination correlation, 
and using the rotational periods from the recent literature, we constrain the stellar inclination for three of the stars. An Li abundance 
analysis yields an independent estimate of the age, but this is inconsistent with the asteroseismically determined age for one of the 
stars. We assess the performance of five grid-based analysis methods and find them all to provide consistent values of the surface 
gravity to ~0.03 dex when both atmospheric and seismic constraints are at hand. The different grid-based analyses all yield fitted val- 
ues of radius and mass to within 2.4cr, and taking the mean of these results reduces it to l.Scr. The absence of a metallicity constraint 
when the average large frequency separation is measured with a precision of 1% biases the fitted radius and mass for the stars with 
non-solar metallicity (metal-rich KIC 11395018 and metal-poor KIC 10273246), while including a metallicity constraint reduces the 
uncertainties in both of these parameters by almost a factor of two. We found that including the average small frequency separation 
improves the determination of the age only for KIC 1 1395018 and KIC 1 1234888, and for the latter this improvement was due to the 
lack of strong atmospheric constraints. 

Key words. Asteroseismology - Stars; solar-type - Stars: fundamental parameters - Stars: atmospheres - Stars: individual: 
KIC 11395018, KIC 10273246, KIC 10920273, KIC 10339342, KIC 11234888 - Stars: oscillations 



1 



1. Introduction 

Stars like the Sun have deep convective envelopes where 
stochastic excitation gives rise to a rich spectrum o f res- 



onant oscillation modes (e.g. Brown & Gillilandl Il994t 



IChristensen-DalsgaardI l2005t lAerts et all l2010l) . The fre- 
quencies of the modes depend on the journey that the 
waves make through the star, so that if the seismic signa- 
tures can be observed, t hey provide v er y accurate probes 
of th e stellar interior (UlrichL 1970| : Leibacher & Stein , 



19711; IChristensen-Dalsgaardl 120071: iMetcalfe et al 



Deheuvels & Michel[. 201 Oat de Meulenaer et al.L 



201C; 



201C; 



Van Grootel et al.L l20ld ICreevev & Bazotl 1201 iTT rhe oscil- 



lations have tiny amplitudes, and with photometry they can 
only be revealed with very long high-preci s ion time series from 



space e.g. 



with CoRoT 



long nig j 
dBadinL 



20031: lAppourchaux et al 



200 8': 'Garcia etal.l, 120091: ' Mosser et al.L 120091: iDeheuvels et al 
2010b: Mathur et al 



The Kepler satellite (iBorucki et al.L 1201 Oh . which was 
launched in early 2009, is providing photometric data 
of outstandin g qual it y on thousan ds of s ta rs, see e.g. 
[Bedding et al.' ('201^, 'Gilliland et al.' ('2010a'), 'Hek ker et al.l 
(i20ia)." .Huberetal.. (,20101 Kallingeret al.. (2010l) . and 



IStello et al.l(l2010h . With an average cadence of 30 minutes, the 
primary objectives of the mission, to search for and characterise 
Earth-like planets, can be reached. For a smaller sample of 
stars (~512), the short cadence of one minute provides the 
time-sampling adequate for detecting oscillation signatures 
present in solar-Uke stars, thus enabling a characterization of 
such planet-hosting stars. 

Due to the low amplitude of the oscillation modes, some 
of the individual frequencies (y/,„ with degree I and ra- 
dial order «) may not be detectable. However, many anal- 
ysis techniques allow one to determine seismic signatures 
in relatively low s ignal-to-noise ratio (S/N) power spectra 
dHuber et all 120091: iMosser & Appourchauxi [20091: [Roxburgh , 



2009; 



201C; 



Campante et all 120101: iHekker et al.L 120101: iKaroff etal 



Mathur et al. . 2010bl) . These signatures are primarily i) 
the average large frequency separation (Av) where Av/,„ = v;.„ - 
v/,„_i , and ii) the frequency corresponding t o the maximum of the 
bell-shaped amplitude spectrum v„,ax (e.g. lHuber et al.ll2010l) . If 
the individual frequencies are available , then (Av) can be d eter- 
mined with somewh at higher precision. IHekker etaLldlonb and 
IVemer et al.l (|201 ll) compare (Av), Vmax, and the uncertainties 
from a variety of established analysis techniques. 

Anticipating the seismic quantities (Av) and Vnvdx to be mea- 
sured for many stars, various pipeline methods based on stellar 
evolution and structure models have been developed to use these 
data to determine stellar properties, such as radius and mass, in 
an efficient and automatic manner In the ideal case, this type 
of grid-based approach provides good estimates of the parame- 
ters for a subsequent detailed seismic study. However, when os- 
cillation frequencies are not available, grid-based methods still 
provide reliable estimates of the mean density, surface gravity, 
radius, and mass. For example, Stello et al. (2009a) compared 
the radius determined from different automatic analyses using 
simulated data a nd find that the r adius can be determined with a 
precision of 3%. lGaietan(l201ll) made a detailed study of grid- 
based methods for asteroseismology using simulated data and 
also solar data. They investigated the errors in the parameters 
and concluded that the surface gravity can be determined with 
practically no systematic bias. Quirion et al. (2010) compared 
their automatic determination of stellar parameters with direct 
measurements of mass and radius for eight bright nearby targets 



(using interferometry and/or binaries), and found agreement to 
within Icr for all stars except one. Determining accurate stellar 
properties, in particular for single stars with V > 7, is an impor- 
tant step towards understanding stellar structure and evolution all 
across the HR diagram, and seismic observations provide possi- 
bly the most accurate way of doing this. This, in t urn, can help i n 
studies such as galactic stellar populations, e.g. iMiglid (1201 ll) : 
.Mo sser et al. (2011) 

Several stars observed by Kepler were selected to be moni- 
tored at the short-cadence rate for the full duration of the mis- 
sion in order to test a nd validate the time series photometry 
(iGilliland et al.Ll2010bh . In this paper we study five of these stars 
which show clear solar-like oscillation signatures. We use the 
global seismic quantities to determine their surface gravity, ra- 
dius, mass, and the age while also assessing the validity of grid- 
based analyses. 

The five stars have the following identities from the 
Kepler Input Catalog (KIC); KIC 11395018, KIC 10273246, 
KIC 10920273, KIC 10339342, and KIC 11234888, hereon re- 
ferred to as CI, C2, C3, C4, and C^l respectively, and their 
characteristics are given in Table [1] The Kepler short-cadence 
Q01234 (quarters 0-4) time series photometry provides the 
global seismic quantities. The atmospheric parameters were de- 
termined from both ground-based spectroscopic data — anal- 
ysed by five different methods — and photometric data (Sect. 2). 
Five grid-based analysis methods based on stellar models are 
presented (Sect. 3) and used to determine log g, the mean den- 
sity, radius, mass, and age of each of the stars (Sect. 4). We also 
test the influence of an additional global seismic constraint, the 
mean small frequency separation ((5vo2) where 6vq„ - vq „ - 
V2,n-i. Possible sources of systematic errors, such as the input 
atmospheric parameters, and using different physics in the mod- 
els are discussed (Sect. 5) and then we assess the performance 
of each grid-based analysis method (Sect. 6). In Sect.|7]we com- 
bine the stellar properties determined in Sect. 4 with published 
results to provide an asteroseismic distance, to constrain the rota- 
tional period, and to estimate the inclination of the rotation axis. 
We also compare the asteroseismic age to that implied from a 
lithium abundance analysis. 

2. Observations 

2.1. Seismic observations 

The Kepler targets CI, C2, C3, C4, and C5 have been observed 
at short cadence for at least eight months (QO-4) since the begin- 
ning of Kepler science operations on May 2, 2009. Observations 
were briefly interrupted by the planned rolls of the spacecraft 
and by three unplanned safe-mode events. The duty cycle over 
these approximately eight months of initial observations was 
above 90%. After 252 days, two CCD chips failed, which af- 
fected the signal for the targets CI and C4. The time series 
were analysed using the raw data provid ed by the Kepler Science 
Operations Center (Jenkins et al.L 2010|). subsequently corre cted 
as described bvlGarcfa et alT(l201 ll) . Campante et al.l (1201 ll) and 
iMathur et al.l (1201 lb presented details on the data calibration, as 
well as an in-depth study of the time series for stars CI, C2, C3, 
and C5 using a variety of documented analysis methods. In this 
paper we use two methods to determine th e global seism i c quan- 
tities: an automatic pipeline package A2Z (iMathur et al.L 1201 Obi) 
and a fit to the individual frequencies. 



' Each of these stars has a pet cat name assigned to it within this 
collaboration. These are Boogie, Mulder, Scully, Cleopatra, and Tigger, 
respectively. 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 
Table 1. Basic data for the five solar-type stars. 



Kepler ID 


Adopted 


Cat name 


RA 


DEC 


Kp 


Time series 


Spectral 




Name 




(hrs min sec) 


(° ' ") 


(mag) 


(days) 


Type* 


KIC 11395018 


CI 


Boogie 


19:09:55 


49:15:04 


10.762 


252.71 


G4-5IV-V 


KIC 10273246 


C2 


Mulder 


19:26:06 


47:21:30 


10.903 


321.68 


F9IV-V 


KIC 10920273 


C3 


Scully 


19:27:46 


48:19:45 


11.926 


321.68 


G1-2V 


KIC 10339342 


C4 


Cleopatra 


19:27:05 


47:24:08 


11.984 


321.68 


F7-8IV-V 


KIC 11234888 


C5 


Tigger 


19:07:00 


48:56:07 


11.926 


252.71 





* Spectral type determined from the ROTFIT method. 



2.1 .1 . Determination of seismic parameters using tine A2Z 
package 

The A2Z pipeline looks for the periodicity of p modes in the 
power density spectrum (PDS) by computing the power spec- 
trum of the power spectrum (PS2). We assume that the highest 
peak in the PS2 corresponds to {Av)/2. We then take a 600 yuHz- 
wide box in the PDS, compute its PS2 and normalise it by the 
standard deviation of the PS2, cr. We repeat this by shifting the 
box by 60 yuHz and for each box, we look for the highest peak 
in the range [<Av)/2 - 10 yuHz, <Av)/2 + 10 juHz]. The boxes 
where the maximum power normalised by cr is above the 95% 
confidence level threshold delimit the region of the p modes, 
[/inin,/max]- The Uncertainty on {Av) is taken as the value of the 
bin around the highest peak in the PS2 computed by taking the 
PDS between f^i„ and f^.^^. 

We estimate the frequency of maximum power, Vmax, by fit- 
ting a Gaussian function to the smoothed PDS between and 
/max- The central frequency of the Gaussian is Vmax, and its un- 
certainty is defined as the smoothing factor, Ix(Av), resulting in 
a precision of between 6% and 8%. 



2.1 .2. Determination of <Av) from individual frequencies 

The frequencies of the oscillation m odes for CI, C2, C3 , and 
C5, have recent l y been published by iMathur et al.l (1201 ih and 
ICampante etaH il20T U). We Ust these frequencies in the ap- 
pendix. (For C4 the S/N in the power spectrum is too low to 
ac curately determine the frequencies.) Following the approach 
bv lWhite et al.l (1201 lb . we performed an unweighted linear least- 
squares fit to the / = frequencies as a function of radial order 
using the available range of frequencies [/lnin,/max]- The fitted 
gradient of the line is (Av), and the uncertainties are derived di- 
rectly from the fit to the data. 

We used (Av) determined from this method together with 
Vmax from the A2Z pipeline as the seismic data for CI, C2, 
C3, and C5. For C4 we used both (Av) and v,nax from the A2Z 
pipeline. Table|2]lists the average se ismic parameters, and t hese 
are in agreement wi th those given in lCampante et al.l (1201 lb and 
iMathur et all (l201lh . 



2.1.3. Mean small frequency separation 

Because the individual frequencies are available, we can read- 
ily calculate the mean small frequency separation (i5vo2), which 
serves as an extra constraint on the stellar parameters. However, 
calculating this quantity from the models implies a calculation of 
the theoretical oscillation frequencies for each model, which is 
not the main purpose of most of the pipelines decribed here, and 
it is generally an observable that is not available for stars with 
low S/N power spectra such as C4. We derived (i5vo2) from the 



Table 2. 

data. 



Mean seismic parameters determined from the Kepler 



KIC ID 


Star 


<Ay> 


^max 


fmm-> fmax 


<(5V02> 






(/uHz) 


O^Hz) 




(/"Hz) 


11395018 


CI 


47.52±0.15 


834±50 


686,972 


4.77±0.23 


10273246 


C2 


48.89±0.09 


838±50 


737,1080 


4.40±0.44 


10920273 


C3 


57.27±0.13 


990±60 


826,1227 


4.76±0.14 


10339342 


C4 


22.50±1.50 


324±25 


219,437 




11234888 


C5 


41.81±0.09 


673±50 


627,837 


2.59±0.40 



individual frequencies within the [/lnin,/max] range and we list 
these values in Table |2] however, we only use it as a constraint 
on the models in Sect. 14.: 



2.1.4. Solar seismic parameters 

We analysed the solar frequencies in the same way as the five 
Kepler stars. We derive (Av) - 135.21 + 0.11 fAiz using the 
frequency range [1000,3900] uH z and t he oscillation f requen- 
cies from iBroomhall et al.l (l2009b . while IMathur et all (1201 Obb 
derived Vniax= 3074.7 + 1.02 /iHz for the Sun. The uncertainty 
in Vm-dx is lower than Ix(Av), and so we artificially increased the 
uncertainty to 145 //Hz (5%) for a more hor nogenous ana l ysis, i n 
line with its expected precision according to lVerner et all (1201 lb . 
We also calculate (5vo2)= 8.56 ± 0.28 //Hz. 

2.2. Spectroscopic Observations 

We observed the targets CI - C4 using the FIBS spectrograph 
on the Nordic Optical Telescope (NOT) telescope located in 
the Observatorio del Roque de los Muchachos on La Palma. 
The targets were observed during July and August 2010 using 
the medium-resolution mode (R - 46,000). Each target was ob- 
served twice to give total exposure times of 46, 46, 60, and 60 
minutes respectively. This resulted in an S/N of ~80, 90, 60, and 
60 in the wavelength region of 6069 - 6076 A. The calibration 
frames were taken using a Th-Ar lamp. The spectra were reduced 
using FIESTOO10- 

The reduced spectra were analysed by severa l groups in- 
dependently using t he following meth ods: SOU (ISousa et al 
2007, .20081), VWA (iBruntt etalfaoiOb 
200( 



20081) . 
, BIA 



ROTFIT (iFrasca et al 



l20Tlh . and NIEM 



dNiemczura & PohibekL l2005b . Here we summarise the main 
procedures for analysing the spectroscopic data, and we refer 
readers to the appendix and the corresponding papers for a more 
detailed description of each method. 



http://www.not.iac.es/instruments/fies/tiestool/FIEStool.html 
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Two general approaches were taken to analyse the atmo- 
spheric spectra, both using the spectral region ~4300-6680 A. 
The first involved measuring the equivalent widths (EW) of 
lines and then imposing excitation and ionization equilibrium 
using a spectroscopic analysis in local thermodynamic equi- 
librium (SOU, BIA). The second approach was based on di- 
rectly comparing the observed spectrum with a library of syn- 
thetic spectra or re ference stars (see, e.g., iKatz et all 119981: 
ISoubiran et al.l[T998h . either using full regions of the spectrum 
(ROTFIT), or regions around specific lines (VWA, NIEM). The 
atomic lin e data were taken fro m the Vienna Atomic Lin e 
Database ('Kuoka et a l.L Il999h and ICastelU & Hubri g ('2004£. 
The MARCS ' (.Gustafsson et all l2008h and ATLAS9 tKurucz , 
model atmospheres were use d, and the syn thetic spectra 
were computed with the SYNTHE (iKuruczl [T993h and MOOG 
(ISne den. 1973) codes. An automatic spectral type classification 
(cf. Table [B is given by ROTFIT. 

The derived atmospheric parameters for CI - C4 for each 
method are given in Table [3] and Fig. [T] shows the fitted 
Tgff and logg for stars CI and C2 (grey and black, respec- 
tively, top panel), and C3 and C4 (grey and black, respectively, 
bottom panel). Each symbol represents the results from one 
spectroscopic analysis: A=SOU, n^ROTFIT, o=VWA, o=BIA, 
x=NIEM. The figures emphasise the correlations between the 
two parameters, especially for CI and C2: a lower T^s is usu- 
ally fitted with a lower logg. The dotted lines represent the 
asteroseismic determination of log g, as explained in Sect. 14.11 
Considering the low S/N of the intermediate resolution spectra, 
we find that there is an overall good agreement (within l-2cr) 
between the methods. However, it is also clear from Table[3]that 
there are some trends corresponding to the method used. For ex- 
ample, the BIA method gives a systematically higher T^s than 
the VWA method, the ROTFIT method generally yields higher 
log g than the other methods, and the metallicity is on average 
higher when determined using EW methods (SOU, BIA) than 
with the line-fitting methods (NIEM, VWA, ROTFIT). A gen- 
eral comparison between the spectroscopic methods based on a 
much larger sample of stars with high resolution spectra is cur- 
rently on-going, and will be reported elsewhere. 



X 
T3 

O 




5400 5600 5800 6000 6200 6400 

Te„ (K) 




X! 
&0 
O 



5400 5600 5800 6000 6200 6400 

Ter, (K) 

Fig. 1. TefF and logg derived from five spectroscopic analyses. 
Each symbol represents the results from one analysis method: 
A=SOU, n^ROTFIT, o=VWA, o^BIA, x=NIEM. The dashed 
lines show the (loggMR) values derived from the seismic data 
alone (see Sect. 14. lb . 



2.3. Photometrically derived atmospheric parameters 



Ground-based SLOAN-griz photometry is available for a large 
number of stars in the field of view. The Kepler Input Catalog 
(KIC) lists the magnitude in the wide Kepler band pass, Kp, 
as well as the griz magnitudes and the stellar parameters (Teff, 
[Fe/H], logg, radius R) derived using these data. However, the 
primary purpose of the KIC was to allow discrimination of 
dwarfs from other classes of stars to aid in the selection of 
planet-hosting candidates. It has become clear since the time se- 
ries data became available t hat the KIC Tpff are not always acc u- 
rate on a star-t o -star b asis ( iMolenda-Zakowicz et al.L l2010allbl : 
iLehmann et al.L l201 ih . T^g were, therefore, re-calculated by 
Pinsonneault & An (201 1 ) using SLOAN photornetry an d the 
YREC models (lAn et al.L l2009t iDemarque et"al1. l2008h . and 
cross-checking the results using the i nfra-red flux method cal - 
ibration based on 2MASS photometry (ICasagrande et al.Ll2010h . 
The Tetf and uncertainties are listed in Table |3] with the heading 
'P/A'. 



3. Seismic methods 

The nearly uninterrupted short-cadence Kepler time series yield 
power spectra that exhibit some signatures of oscillations. Even 
from low S/N power spectra, the global seismic quantities (Av) 
and Vmax can easily be determined without having to measure 
individual frequencies. In preparation for the hundreds of stars 
in the Kepler data where these seismic parameters are read- 
ily available, several pipeline codes have been developed to 
use stellar models to infer the surface gravity log g, mean den- 
sity {p}, radius R, mass M, and age t of stars from the global 
seismic quantities supplemented with atmospheric parameters, 
such as Teff. In this analysis we used the followi ng methods 
which are discussed belo w: CESAM2k/Mumbai (Mazumdail 
l2005h . Yal e-Birmingham (iBasu etakL |201(A lOaietakL 1201111 



RADIUS (IStello et al.L l2009aL iMetcalfe et all 120 ioT SEEK 



(lOuirion et al. . 20IOI) . and RadExlO (Creevey in prep.). We used 



http://www.user.oat.ts.astro.it/castelli/grids.html 



five different analysis methods in order to test the validity of our 
results, and to assess the performance of each code for producing 
reliable stellar parameters. 
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Table 3. Atmospheric parameters determined by various spectroscopic analyses of the same NOT spectra, and from photometric 
analysis of the KIC photometry ('P/A'). 





T 

cfF 


logs 


Lre/HJ 




V sin f 




(K) 


(dex) 


(dex) 


(kms"') 


(kms-') 


CI 












sou 


5717+68 


3.96+0.11 


+0.35+0.05 


1.30+0.03 




ROTFIT 


5445±85 


3.84+0.12 


+0.13+0.07 




1.1+0.8 


VWA 


5580±79 


3.81+0.12 


+0.19+0.06 


1.40+0.13 




BIA 


5650+60 


4.10+0.10 


+0.36+0.10 


1.20+0.20 




NIEM 


5700+100 


4.10+0.20 


+0.13+0.10 


0.70+0.40 






DooU+D / 










C2 












SOU 


6165+77 


4.01+0.11 


-0.04+0.06 


1.48+0.05 




ROTFIT 


5933+205 


4.07+0.10 


-0.21+0.08 




3.2+1.5 


VWA 


6050+100 


3.80+0.11 


-0.18+0.04 


1.50+0.10 




BIA 


6200+60 


4.00+0.20 


-0.04+0.07 


1.50+0.20 




NIEM 


6200+100 


3.90+0.20 


-0.18+0.05 


0.50+0.40 




r/A 


OJOU+ lb 










C3 












sou 


5770+75 


4.08+0.11 


+0.04+0.05 


2.11+0.08 




ROTFIT 


5710+75 


4.15+0.08 


-0.02+0.07 




1.5+2.2 


VWA 


5790+74 


4.10+0.10 


-0.04+0.10 


1.15+0.10 




BIA 


5800+60 


4.10+0.20 


+0.03+0.07 


1.20+0.20 




NIEM 


6000+100 


3.80+0.20 


-0.03+0.08 


1.00+0.40 




P/A 


jo<5U+Jo 










C4 












sou 


6217+82 


3.84+0.11 


-0.11+0.04 


1.60+0.20 




ROTFIT 


6045+125 


4.03+0.10 


-0.23+0.08 




4.0+2.8 


VWA 


6180+100 


3.65+0.10 


-0.15+0.10 


1.75+0.10 




BIA 


6200+100 


3.70+0.20 


-0.06+0.08 


1.60+0.20 




NIEM 


6200+100 


3.70+0.20 


-0.17+0.06 


0.50+0.40 




P/A 


6280+63 










C5 












P/A 


6240+60 











3.7. CESAM2k/Mumbai 

The analysis usin g the CESAM2k stellar evolution code 
(iMorel & Lebretoni l2008h is based on the comparison of 
both seismic and non-seismic observations (^={(Av),Vtnax, 
Tgf{,logg,[Fe/H],{6vo2}}) with those calculated from a grid of 
stellar evolution models. This version of CESAM2k uses the 
OPAL equation o f state ( Rogers & Navfonovl l2002h and the 
OPAL opacities dlglesias & Rogers[~ 19961) supplemented b y 
the low temperature opacities of Alexander & Ferguson ( 1994h . 
The solar mixture is given by iGreye sse & Sauvaf 71998^" and 
the N ACRE nuclear reaction rates (lAngulo. Arnould. & Ravei 
[1991 are used . Convection is described by the standard mixing 
length theory jBohm-Vitensel Il958h . The models also include 
microscopic diffus ion of helium and h e avy el ements, following 
the prescription of iProffitt & MichaudI (|1991|) for masses < 1.3 
Mo. 

The grid of models used in this analysis spans the mass range 
of 0.80 Mo to 1.70 Mo, in steps of 0.02 Mo. The initial metal- 
licities of the models Zi range from 0.005 to 0.030 in steps of 
0.005, and for each value of Zj, five different combinations of 
(Xi, Yi) are used, where Y, is the initial helium fraction. Three 
mixing length parameters a = 1.8, 1.9,2.0 are considered, and 
for stars with convective cores, we used convective overshoot 
to the extent of x Hp (Hp - pressure scale height), with 



ffov = {0.00, 0.05, 0.10, 0.15, 0.20, 0.25). The stellar evolutionary 
tracks start from the zero age main sequence (ZAMS) and con- 
tinue until log(reff) ~ 3.715, which corresponds to on or near 
the red giant branch. For each model during the evolution, os- 
cillation frequencies for low degree modes {£ - 0, 1 , 2, 3) are 
computed under the adiabatic approximation using the Aar hus 
pulsation package, ADIPLS (Christensen-Dalsgaardi,l2008bl) . 

For every star, the average large frequency separation was 
computed between the same frequency limits as used in the ob- 
served data for the corresponding star and run. The average was 
calculated by integrating the individual frequency separation as 
a function of frequency between the given limits and dividing 
by the frequency range. The frequency of maximum amplitude, 
Vm-dx, has been calculated using the scaling relation provided by 
iKjeldsen & Bed ding (1995.) : 



M/Mo 



(R/Rey- V7W5777K 



3,050yuHz 



(1) 



lows: 



We constructed a reduced using the available data as fol- 



X 



"Diu- , ^2 



z 



O" /,obs 



(2) 
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where the sum is over all the available constraints and xj{ 
takes into account the number of constraints. The subscripts 
"obs" and "mod" refer to the observed and model values, respec- 
tively. The quoted uncertainty in the ith observed data is denoted 
by cr,- obs. We sought the minima of^^ to provide estimates of the 
stellar parameters. The quoted values of the parameters and their 
uncertainties are the midpoints and half the span of the ranges of 
parameters for which < 1. 

3.2. Yale-Birmingham 



3.3. RADIUS 



The Yale-Birmi ngham (iBasu et al. I I20T0I) code as described by 
iGai et al.] ( 1201 lb is a grid-based method for determining a star's 
mass, radius, and ag^ The grid of models was constructed 
using the Yale Rotation and Evolution Code (iDemarque et al.L 
l2008h in its non-rotating configuration. The input physics in- 
cludes the OPAL equation of state tables, the OPAL temper- 
ature opacities supplemented with low tem perature opacities 
from [Ferguson. Alexander. & AllardI ( l2005h . and the NACRE 
nuclear reaction rates. All models include gravitational set- 
tling of helium and heavy ele ments using the formulation of 
iThoul. Bahcall. & Loebid 19941) . We use the Eddington ^(t) rela- 
tion (here r means optical depth), and the adopted mixing length 
parameter is a - 1 .826. An overshoot of q-qv = 0.2 was assumed 
for models with convective cores. The grid consists of 820,000 
individual models with masses ranging from 0.8 to 3.0 Mq in 
steps of O.2M0. These models have [Fe/H] ranging from -1-O.6 
to -0.6 dex in steps of 0.05 dex. We assume that [Fe/H] = 
correspon ds to the solar abundance (Zq/Xq - 0.023) as deter- 
mined by IGrevesse & Sauvall ( Il998h . The model value for v^ax 
is calculated using Eq. ([Tj), and (Av) is determined using 



(Av) 
<Ay)o 




(3) 



where <Av>0 = 134.9 yuHz (see lKieldsen & Bedding|ll995h . 

The Yale-Birmingham pipeline finds the maximum likeli- 
hood of a set of input parameters calculated with respect to the 
grid of models. The estimate of the parameter is obtained by 
taking an average of the points that have the highest likelihood. 
We average all points with likelihoods over 95% of the maxi- 
mum value of the likelihood functions. For a given observational 
(central) input parameter set, the first key step in the method is 
generating 10,000 input parameter sets by adding different ran- 
dom realizations of Gaussian noise to the actual (central) obser- 
vational input parameter set. The distribution of any parameter, 
say radius, is obtained from the central parameter set and the 
10,000 perturbed parameter sets form the distribution function. 
The final estimate of the parameter is the median of the distri- 
bution. We used Icr limits from the median as a measure of the 
uncertainties. 

The likelihood function is formally defined as 



n 



1 



V2^, 



,7TCri 



(4) 



where is given in Eq.|2] and q = {Teff, [Fe/H], Av, Vmax). From 
the form of the likelihood function in Eq.|4|it is apparent that we 
can easily include more inputs, or drop some inputs depending 
on the availability of data. 



The RADIUS pipehne dStello et all l2009al) is based on a large 
grid of ASTEC models (fchristensen-Dalsgaarc , 1200 8a') using 
the EFF equation of state ( Eggleton et al.L 11973 ). We u s ed the 
opacity tables of iRogers & IglesiasI ( 1995|) and Kuruczl (Il99ll) 



(for r < 10 K), with solar mixture of lGrevesse & Noel; 



Rotation, overshooting, and diffusion were not included. The 
grid was created with fixed values of the mixing-length parame- 
ter, a — 1.8, and the initial hydrogen abundance of X\ - 0.7. The 
resolution in log(Z/X) was 0.1 dex between 0.001 < Z < 0.055, 
and the resolution in mass was 0.01 M© from 0.5 to 4.0 Mq. The 
evolution begun at the ZAMS and continued to the tip of the red 
giant branch. To convert between the model val ues of Z and the 
observed [Fe/H], the pipeline used Zq = 0.0188 (ICoxLl2000l) . 

Each output parameter was determined by selecting the set 
of models that were within +3cr of the observed input data. We 
pinpointed a single best-fitting model using a formalism and 
the Icr uncertainty is estimated as 1/6 of the maximum range of 
the valu es of the selected m odels. The pipeline as described in 
detail bv lStello et all ( l2009ab had some slight modifications; for 
example, the large frequency separation was derived by scaling 
the solar value (see Eq. O instead of calculating it directly from 
the model frequencies. 



3.4. SEEK 



In the absence of a metallicity measurement the determination of 



In tne absence or a metallicity mea 
the age is hindered, see lGaietalJ ( l201lh 



The SEEK procedure dOuirion et a ll 120 lot) also makes use 
of a large grid of stellar models computed with the ASTEC 
code. This version of ASTEC uses th e OPAL equation of 
state (iRogers. Swenson. & IglesiasLll996l) al ong with the OPAL 
plus Fer guson & Alexander opacity tables (llglesias & RogersL 
1 1996; Alexander & Fergusonl 1 19941). the element to element ra- 
tios in the metallic mixture of lGrevesse & Sauva I (119981) . and 
co nvection is treated with the mixing-length formulation of 
Bo hm- Vitensel d 1 95 8l) ; the mixing length to pressure scale height 
ratio a, characterizing the convective efficacy, is treated as a vari- 
able parameter in the SEEK fits. Neither diffusion nor overshoot- 
ing is included. Oscillati on frequencies for each model ar e calcu- 
lated using the ADIPLS (IChristensen-Dalsgaardl l2008bl) code. 

Two subgrids were created: the first subgrid comprises tracks 
with all combinations of Z = [0.005, 0.01, 0.015, 0.02, 0.025, 
0.03], Xi = [0.68, 0.70, 0.72, 0.74], and a = [0.8, 1.8, 2.8] 
while the second subset has Z = [0.0075, 0.0125, 0.0175, 0.0225, 
0.0275], Xi = [0.69, 0.71, 0.73], a = [1.3, 2.3]. Every subset is 
composed of 73 tracks spanning from 0.6 to 1.8 M© in steps of 
0.02 M0 and from 1.8 to 3.0 M© in steps of 0.1, and each track 
was evolved until just after the base of the giant branch or t = 
15x10^ yrs. The metallicity combinations correspond to -0.61< 
[Fe/H] < 0.20. 

SEEK compares an observed star with every model of the 
grid and makes a probabilistic assessment of the stellar param- 
eters, with the help of Bayesian statistics. Its aim is to draw the 
contour of good solutions which is located axound x\. The pri- 
ors used in that assessment are flat for the age, the metallicity, the 
initial helium ratio, and the mixing length parameter. The only 
non- flat prior is relate d to the initial mass function and makes use 
of the|Chabrii3d2Q01]) IMF model, where ^(M) = 0.019M", with 
n = -1.55 for M < I.OMq and n = -2.70 when M > I.OMq. 
The details of the SEEK procedure, including the choice of pri- 
ors, and an introduc tion to Bayesian statistics can be found in 
lOuirion etaDd2010l) . 
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3.5. RadExW 

RadExlO is a grid-based approach to determining the radius, 
mass, and age using some or all of the following as input 
{(Av),Vn,ax,reff,log^,[Fe/H]). It is based on the ASTEC code and 
uses the EPF equation of state of lEggleton et al.l (Il973l) without 
Coulo mb corrections, the OPAL opacities "IS esias & RogersL 
1 19961) supplemented by Kurucz opacities at low temper- 
atures, and solar mixture from iGrevesse & Noels' ('19931). 
The n uclear reaction rates came from Bahcall & Pinsonneauij 
( 19921). convection is described by the mixing-length theory of 
lBohm-Vitens3 (119581) . convective core overshooting is included 
with aov set to 0.25, and diffusion effects are ignored. 

The grid considers models with masses from 0.75 - 2.0 
in steps of 0.05 M©, ages fromZAMS to subgiant, Zi spans 0.007 
- 0.027 in steps of ~ 0.003, while Xi is set to 0.70: this cor- 
responds to Yi = 0.263 - 0.283. The mixing length parameter 
a = 2.0 is used, which was obtained by calibrating the solar data. 
To obtain the stellar properties of mass, radius, and age, we per- 
turb the observations using a random Gaussian distribution, and 
compare the perturbed observations to the model observables to 
select an optimal model. The R, M, t, and their uncertainties are 
defined as the mean value of the fitted parameter from 10,000 
realizations, with the standard deviations defining the Icr uncer- 
tainties. 



4. Stellar properties 

4. 1 . Constraining log g witli seismic data 

The quantity {Ay) is proportional to the mean density of the star 
(see eq. [S)) and Vmax also scales with R and M (see eq. IH). 
By making an assumption that these stars have roughly solar T^ff 
with a large uncertainty, the seismic data alone should give a 
robust estimate of log g by using the scaling relations for v^ax 
and (Av) and solving for R and M. The first two data columns in 
Table ID show \ogg and its uncertainty, denoted by the subscript 
'v', when we use the seismic data and a T^ff estimate of 6000 K 
+ 500. 

The grid-based methods described in Sect. [3] also used {Av) 
and Vtnax from Table [2] as the only input observational data to 
their codes to obtain a model asteroseismic value of log g. In this 
case we restricted the model Teff to less than 8000 K. Table |4] 
shows the mean value of log^, (loggMR), obtained by combin- 
ing the results from the five methods. We indicate that these are 
model-determined values by the subscript 'MR' . As can be seen, 
there is very good agreement between the model-independent 
and model-dependent values of log g. 

Figure |2] shows the difference between the asteroseismic 
\ogg value returned by each method, logg,, and (log^MR)- 
From left to right on the x-axis we show the results 
from SEEK, RADIUS, RadEXlO, Yale-Birmingham, and 
CESAM2k/Mumbai (note the abbreviated labelling on the x- 
axis). For each seismic method, we show the differences in \ogg 
from left to right for stars CI, C2, C3, C4, C5, and the Sun. For 
the CESAM2k/Mumbai method 'CMum', we show results for 
CI, C2, and C3 only, and these are labelled accordingly to avoid 
confusion. 

The fitted log g values for each star and each method differ 
by ~0.05 dex, excluding C4, with the largest differences seen 
primarily between the RADIUS and Yale-Birmingham meth- 
ods, but still matching within Icr. For C4 SEEK and RadExlO 
give the largest dispersions, this difference reaching 0.13 dex. 
However, C4 is a very evolved MS star (logg~3.50 dex) and 



Table 4. Surface gravity obtained from (Av) and Vmax using scal- 
ing relations and stellar models, denoted by subscripts 'v' and 
'MR', respectively. 



Star 




(Ty 


(loggMR)" 


O-sys 




" Stat 




(dex) 


(dex) 


(dex) 


(dex) 


(dex) 




CI 


3.88 


0.02 


3.88 


0.05 


3.87 


0.07 


C2 


3.88 


0.03 


3.89 


0.04 


3.89 


0.06 


C3 


3.95 


0.03 


3.98 


0.05 


3.97 


0.06 


C4 


3.47 


0.04 


3.49" 


0.13" 


3.42 


0.06 


C5 


3.78 


0.04 


3.82 


0.05 


3.80 


0.07 


Sun 






4.42 


0.02 


4.43 


0.02 



"(log^MR) is the mean value of logg provided by all of the model re- 
sults. 

'log gMR and (Tsta, are the SEEK values. 

Discarding the value from RadExlO yields (logg) = 3.46 and a^y^ = 
0.07 dex. 



the results for RadExlO are biased, since this grid concentrates 
mostly in the MS and just beyond. If we ignore the result for 
RadExlO, we find a maximum difference between the results of 
0.07 dex. The differences of 0.05 and 0.07 dex are still in much 
better agreement than the spectroscopic methods. 

We also see that the statistical errors given by each method 
vary between 0.01 - 0.05 dex. However, all of the results fall to 
within 2o" of the mean value. Each reported uncertainty is not un- 
derestimated, but some grids take into account some other vari- 
ables, e.g. different mixing-length parameters, that will increase 
the reported value. For example, the uncertainties reported by 
Yale-Birmingham and RADIUS are indeed correct, however, be- 
cause their grids are based on different physics and sets of pa- 
rameters, they obtain results that differ by more than Icr. Gai 
et al. (2010) did a systematic study of the uncertainties in stel- 
lar parameters determined by grid-based methods, including the 
Yale-Birmingham grid. Their study included a detailed investi- 
gation of the errors in logg, and they found typical statistical 
uncertainties of 0.014 dex when Tgff is included, in agreement 
with those reported here. They also tested the systematic errors 
by using different grids and different parameters (e.g. a differ- 
ent mixing-length parameter) and found that for stars similar to 
those used in this paper i.e. with similar (Av), the half-width at 
half-maximum of the distribution of systematic errors is around 
3% in logg (see their Fig. 19, panel b). For a logg value of 
around 3.8, this would imply cr < 0.1 dex, in agreement with 
the differences in logg reported above. 

As another example, the CES AM2k/Mumbai method reports 
results that take into account different values of the convective 
core overshoot parameter in the models. The inclusion of more 
parameters/physics will yield a more conservative error. 

While distinguishing between the statistical uncertainties 
and the systematic errors is not always clear, we can be confident 
that combining the results from different grids provides a reliable 
determination of log g while the dispersion among the results is 
representative of a typical systematic error found by using dif- 
ferent sets of physics and parameters. In order to provide accu- 
rate results with a conservative error we adopt the results from 
SEEK, whose uncertainty is larger than the dispersion among the 
fitted results. In Table|4] last three columns, we give an estimate 
of the systematic error, cr^y^ - max{| loggsEEK - logg,|), and the 
log g values and uncertainties provided by SEEK. 
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T3 



60 

o 



A 

60 



60 

o 




SEEK RADIUS RadEX YaleBI CMum 

Fig. 2. Comparison of each seismic method's asteroseismic val- 
ues of log g to the mean value (log ^mr) of the four or five meth- 
ods. Each method is abbreviated and labelled on the x-axis, and 
for each method the results are shown from left to right for CI, 
C2, C3, C4, C5, and the Sun. The CMum (CESAM2k/Mumbai) 
method shows results for CI, C2, and C3 only. 



4.2. Radius and mass 

The atmospheric parameters Teff and [Fe/H] are needed to derive 
the radius and mass of the star Using the log g values obtained 
in Sect. 14.11 we selected one set of spectroscopic constraints as 
the optimal atmospheric parameters, and to minimise the effect 
of the correlation of spectroscopically derived parameters. We 
chose the set whose log g values matched closest globally to the 
asteroseismically determined ones. By inspecting Tables |3]and|4] 
we found that the VWA method has the overall closest results in 
log g. We therefore combined these spectroscopic data with (Av) 
and Vn,ax from Table |2] and used these as the observational input 
data for the seismic analysis. For C5 we used the photometric 

In Fig. |3]we show the deviation of the fitted radius of each 
method R/ from the mean value {R ) in units of %, with the mean 
radius in units of given in Table |5] The representation of the 
results is the same as in Fig. |2] i.e. from left to right for each 
method, we show CI, C2, C3, C4, C5, and the Sun. Most of 
the results are in agreement with the mean value at a level of 
Icr, and the uncertainties vary between approximately 1 and 4 % 
(see Sect. 14.11 ). To be consistent, we adopt the SEEK method 
as the reference one. This choice is also justified by the fol- 
lowing reasons: i) SEEK covers the largest parameter space in 
terms of mass, age, metallicity, and mixing-length parameter, ii) 
it has also been tested with direct measurements of mass and 
radius of nearby stars, iii) and it determines a best model param- 
eter for each property (e.g. luminosity, initial metal mass frac- 
tion, and uncertainties), which allows us to make further infer- 
ences as well as investigate the systematics. In addition, it does 
not use the scaling relations, which may introduce a systematic 
bias on the order of 1% in the stellar parameters e.g. lStello et alj 
(l2009bl) . We used the results from the other pipelines as a test of 
the systematic errors. 

In Table|5]we list the R values found using SEEK, the uncer- 
tainties cr, the mean value of R obtained by combining the results 
from all of the grids {R ), osys - max{|/?sEEK - given in units 
of R©, and osys normalised by o",- i.e. max{|(/?sEEK - ^/O/o"/!)- 
Using SEEK with its reference set of physics, we determine the 
radius of each star with a typical statistical precision of 3%. If we 
use the results from the other pipelines as a measure of the sys- 
tematic error, then we report an accuracy in the radius of between 
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Fig. 3. Comparison of each seismic method's determination of R 
with the mean value {R ) using both spectroscopic and seismic 
constraints. See caption of Fig.|2]for details. 



Table 5. Radius determination from the SEEK method (first half 
of table) and combining all methods (second half of table) using 
spectroscopic and asteroseismic data. 



Star 


R 


cr 


cr 


(R) 


CTsys 








(Rg) 


(Rg) 


(%) 


(Rg) 


(Rg) 


(%) 




CI 


2.23 


0.04 


2 


2.21 


0.06 


3 


2.3 


C2 


2.11 


0.05 


2 


2.19" 


0.22° 


11 


2.4 


C3 


1.90 


0.05 


3 


1.88 


0.07 


4 


1.2 


C4 


3.81 


0.20 


5 


3.95* 


0.38 


10 


2.2 


C5 


2.44 


0.14 


6 


2.47 


0.12 


5 


1.0 


Sun 


0.99 


0.03 


3 


1.00 


0.01 


1 


1.0 



" discarding the RADIUS result yields a value of 2.15 Rg ± 0.09. 
* discarding the value from RadExlO yields (R) = 4.0. 



3 and 5% for CI, C3, C5, and the Sun. The RADIUS method re- 
ports a radius that differs by 11 % from the radius of the SEEK 
method for C2, but this difference corresponds to only 2.4crfrom 
the fitted R. Without this value oSys reduces to 5%. However, 
we have no justification to remove this value or discard the pos- 
sibility that this is closer to the correct value. We remove the 
RadExlO value for C4, but this results in an insignificant change 
(< 1%). The method responsible for the largest deviation from 
the SEEK results for C4 is the Yale-Birmingham method, al- 
though they agree just above their Itr level. In the final column 
of the table we show the maximum deviation of the fitted radius 
from the SEEK radius but in terms of the uncertainty given by 
each pipeline method cr,-. Here we see that all of the results are 
consistent with SEEK to 2.4cr. 

Table |6]reports the values of the mass determined by SEEK, 
and the other seismic methods, using the same format as Table|5] 
The statistical uncertainties are of the order of 8%, and the sys- 
tematic errors are typically 5-12%, with each grid-based method 
reporting the same fitted mass to within 2.4cr. The fitted mass 
is highly correlated with the fitted radius, and thus most of the 
same trends are found among each pipeline method for both ra- 
dius and mass. A much larger uncertainty in mass is reported for 
C5, and this is due to the lack of a metallicity constraint. 
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Table 6. Mass determination from the SEEK method (first half 
of table) and combining all methods (second half of table) using 
spectroscopic and asteroseismic data. 



Table 7. Surface gravity and stellar mean density obtained from 
scaling relations (V) and stellar models ('MR') using both seis- 
mic and atmospheric constraints. 



Star 


M 


a 


a 


<M> 


CTsys 


0"sys 




Star 




log gMB. 


ip). 


<P>MR 




(Mo) 


(Mo) 


(%) 


(Mo) 


(Mo) 


(%) 






(dex) 


(dex) 


(kgm-3) 


(kg mr^) 


CI 


1.37 


0.11 


8 


1.35 


0.09 


7 


2.0 


CI 


3.86 ± 0.03 


3.88 ± 0.02 


175 ± 2 


174 ±4 


C2 


1.26 


0.10 


8 


1.37" 


0.43" 


34 


2.4 


C2 


3.88 ± 0.03 


3.88 ± 0.02 


185 ± 1 


189 ±2 


C3 


1.25 


0.13 


11 


1.23 


0.06 


5 


1.0 


C3 


3.94 ± 0.03 


3.97 ± 0.03 


254 ± 1 


257 ±5 


C4 


1.79 


0.12 


7 


1.84* 


0.09* 


5 


1.5" 


C4 


3.47 ± 0.03 


3.52 + 0.04 


39± 16 


46 ±4 


C5 


1.44 


0.26 


18 


1.50 


0.18 


12 


1.6 


C5 


3.79 ± 0.04 


3.82 ± 0.03 


135 + 1 


140 + 2 


Sun 


0.97 


0.06 


7 


1.01 


0.08 


9 


1.0 


Sun 




4.43 ± 0.01 




1400 ± 18 



" discarding the RADIUS result yields a value of 1.32 Ro 
* discarding the value from RadExlO. 



:0.14. 



4.3. log g and (p) using combined seismic and spectroscopic 
data 

Combining {Av) and Vmax with the VWA spectroscopic values 
of Teff and the photometric value for C5, we calculated model- 
independent values of log g just as explained in Sect. 14.11 We 
also calculated the mean density of the star using (Av). Table Q 
lists these properties, again denoted by ' v' to indicate that they 
are obtained directly from the data. We also give the model- 
dependent values of log^MR and {Av )mr as reported by SEEK. 

Both the log g and (p) values are in good agreement when 
derived directly from the data and when using stellar models. 
For \ogg there is agreement to 1.5cr, and for (p) the agreement 
is within 2.5cr. These values represent a relative precision of 2% 
for (p), with the exception of C4. In this case, the model (p) 
provides a more precise value than when using the seismic data 
alone, because the inclusion of atmospheric data helps to nar- 
row down the range of possible values, when the uncertainty in 
(Av) is large. The calculated solar value of 1,400 kg m"-' is in 
excellent agreement with the true solar value (1,408 kg m"^). 

Comparing log gy from Tables |4] and Q we find values that 
differ by at most 0.02 dex (< Icr). This implies that by making 
a reasonable assumption about T^ff, log g can be well estimated 
using only (Av) and Vmax, or Vmax alone. 

Again, inspecting \ogg from Tables |4]and|7] but this time us- 
ing the SEEK model-dependent values ('log^MR') we also find 
that, apart from C4, the derived log g values are consistent with 
and without the atmospheric constraints, but the uncertainties in 
the latter reduce by a factor of 2 - 3. For C4 there is a difference 
of 0. 1 dex between the two values, a change of nearly 2cr. As 
noted before this discrepancy is due to the relatively larger un- 
certainty in (Av) when no atmospheric constraints are available. 

Comparing SEEK to the other pipeline results for \ogg, we 
find o-sys of 0.02, 0.05, 0.01, 0.04, 0.02, and 0.02 dex, respec- 
tively. For the Kepler stars, these differences are smaller than 
o"SEEK except for C2. The agreement between the \ogg val- 
ues from the five seismic methods supports the results that we 
present, and we can be confident that any of the pipeline results 
using both seismic and atmospheric paramete rs can determine 
logg with a minimal systematic bias, just as iGai et"al ] (1201 Ih 
showed. 

We also find that the VWA spectroscopic values of log^ 
agree with the model-dependent ones to within Icr for CI and 
C2, and within 1.5cr for C3 and C4. 



4.4. Age 

Determining the ages of the stars is a more complicated task 
since the value of the fitted age depends on the fitted mass and 
the description of the model. Figure|4]shows the fitted mass ver- 
sus fitted age for stars CI and C3 using five and four seismic 
methods, respectively (the CESAM2k/Mumbai method did not 
fit these data for C3 and C4). For both stars a coiTelation be- 
tween the two parameters can be seen; a lower fitted mass will 
be matched to a higher age and vice versa. However, the uncer- 
tainties from the SEEK method (~ 1 Gyr for a mass of 1.37 Mq 
and ~ 2 Gyr for a 1.25 M© star) do capture the expected correla- 
tion with mass. 

It is known that one of the physical ingredients to the models 
that has a large effect on the age of the star is the convective core 
overshoot parameter (see Sect. l5.2l below) where fuel is replen- 
ished by mixing processes that in effect extends the lifetime of 
the star. The CESAM2k/Mumbai method includes several val- 
ues for this parameter among its grid of models, and this method 
fits an age of 5.3 Gyr for CI (the SEEK method fits 3.9 Gyr). 
However, the fitted mass of the star is also lower than the SEEK 
one, and if we consider the uncertainty arising from the corre- 
lation with the age, then 5.3 Gyr is not outside of the expected 
range. For C2 the fitted mass for CES AM2k/Mumbai is 1 .20 Mq 
(SEEK = 1.26 Mo), and the fitted age is 4.4 Gyr (SEEK = 3.7 + 
0.7 Gyr). The difference between the fitted ages is Icr, again not 
showing any difference due to the adopted values of the convec- 
tive core overshoot parameter 

The final value of the age depends on the adopted model 
(stellar properties and physical description), but what seems to 
emerge from these data is that all of the stars are approach- 
ing or have reached the end of the hydrogen burning phase. In 
Table[8]we summarise the stellar properties obtained by adopting 
(Av) and v^ax from Table |2] the VWA spectroscopic constraints 
from Table |8] and the photometric T^ff for C5, while using the 
SEEK method (left columns) and combining the results from 
four (C3, C4, C5, Sun) and five (CI, C2) grid-based methods 
(right columns). We see that the mean values of the combined 
seismic methods results in stellar properties consistent within 
1.5crof the SEEK results. 



4.5. Including {6vo2) in the seismic analysis 

With the individual oscillation frequencies available, the value 
of {6vq2} can be easily calculated. This observable can be partic- 
ularly useful for determining the stellar age because it is deter- 
mined by the sound-speed gradient in the core of the star, and as 
hydrogen is burned, a positive sound-speed gradient builds up. 
Three of the pipeline codes presented here do not allow one to 
use this value because it relies on the calculation of oscillation 
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Table 8. Stellar properties obtained using (Av), Vm-dx, and the atmospheric constraints from VWA for CI - C4 and P/A for C5, from 
SEEK (first half of table) and by combining the results from all of the pipelines (second half of table). 



Star 


logs 


R 


M 


r 


<logg) 


(R) 


<M) 




(dex) 


(Ro) 


(Mo) 


(Gyr) 


(dex) 


(Ro) 


(Mo) 


CI 


3.88± 0.02 


2.23± 0.04 


1.37+ 0.11 


3.9+ 1.4 


3.87 


2.21 


1.35 


C2 


3.88± 0.02 


2.11+ 0.05 


1.26+ 0.10 


3.7+ 0.7 


3.89 


2.19 


1.37 


C3 


3.97± 0.03 


1.90± 0.05 


1.25+ 0.13 


4.5+ 1.8 


3.97 


1.88 


1.23 


C4 


3.52± 0.04 


3.81+0.19 


1.79+ 0.12 


1.1+0.2 


3.52 


3.95 


1.87 


C5 


3.82± 0.03 


2.44+ 0.14 


1.44+ 0.26 


2.6+ 0.9 


3.82 


2.47 


1.50 


Sun 


4.43+ 0.01 


0.99+ 0.03 


0.97+ 0.06 


9.2+ 3.8 


4.43 


0.99 


0.97 




Table 9. Comparison between fitted stellar properties using 
SEEK without (top lines) and with (bottom lines) {SV02} as a 
seismic constraint. 



1.15 
Fig. 4. Fitted aj 



1.20 



1.40 



1.25 1.30 1.35 
Mass (Ms) 

;e versus mass for stars CI (o) and C3 (□) 



1.45 



frequencies for every model in the grid, unlike (Av) which can 
be derived from a scaling relation (c.f. Eql3]l. We included {SV02} 
as an extra constraint in SEEK to study the possible changes in 
the values and uncertainties of the stellar properties. 

In Table|9]we list log g, R, M, t, and the uncertainties without 
(first line for each star) and with (second line for each star) {6vq2} 
as an observational constraint. The first line for each star shows 
the same values as those presented in Table [8] and we repeat it 
here to make the comparison easier for the reader For CI we 
obtained a more precise determination of the age. For C2 and 
C3, no improvements in the parameters were found. For C5, all 
of the parameter uncertainties reduce by a factor of three, and the 
values of the parameters change by Icr. The usefulness of {SV02} 
here can most probably be explained by the lack of a metallicity 
constraint on the models. For the Sun we also found that the age 
improved from an estimate of 9.2 Gyr to 4.9 Gyr, very close to 
the accepted value. This last result can be explained by the fact 
that the stellar 'observables' change on a much slower scale than 
for a more evolved star and hence provide weaker constraints. 



5. Sources of systematic errors 

5.1. Using different spectroscopic constraints 

We adopted the VWA spectroscopic constraints because their 
log g values were in best agreement with the asteroseismic ones 
given in Table |4] and later confirmed in Sect. 14.31 However, 
the VWA parameters may not be the optimal ones, and so we 
should investigate possible systematic errors arising from using 
the other atmospheric parameters from Table |3] We repeated the 
analysis using {Av) and v^ax with each set of atmospheric con- 
straints to derive radius and mass. In Table [TOl we give the six 
sets of observational data for star CI, and we label these SI - 



Star 


^ogg 
(dex) 




R 

(R0) 




M 

(Mo) 


T 

(Gyr) 


CI 


3.88+ 
3.87+ 


0.02 
0.02 


2.23+ 0.04 
2.22+ 0.05 


1.37+ 0.11 
1.34+ 0.11 


3.9+ 1.4 
4.5+ 0.5 


C2 


3.88+ 
3.88+ 


0.02 
0.02 


2.11 + 
2.11 + 


0.05 
0.05 


1.26+ 0.10 
1.25+ 0.10 


3.7+ 0.7 
3.7+ 0.6 


C3 


3.97+ 
3.97+ 


0.03 
0.02 


1.90+ 
1.90+ 


0.05 
0.06 


1.25+ 0.13 
1.23+ 0.11 


4.5+ 1.8 
5.0+ 1.9 


C5 


3.82+ 
3.85+ 


0.03 
0.01 


2.44+ 0.14 
2.56+ 0.05 


1.44+ 0.26 
1.71+ 0.09 


2.6+ 0.9 
1.6+ 0.2 


Sun 


4.43+ 
4.43+ 


0.01 
0.01 


0.99+ 
1.01+ 


0.03 
0.03 


0.97+ 0.06 
1.01+ 0.04 


9.2+ 3.8 
4.9+ 0.5 



S6. We also give the fitted radius and mass with uncertainties for 
each of these sets using the SEEK method. 

Inspecting sets SI - S5 we see that the largest differences 
between the fitted properties is 0.05 and 0. 12 Mq{~Io-), using 
sets SI and S2. These differences can be attributed to both the 
input Teff and [Fe/H] each contributing about equal amounts. At 
the level of precision of the spectroscopic log g, this value has 
very little or no role to play in determining the radius and the 
mass, due to the tight constraints from the seismic data. 

The results using S6, however, change by more than Icr. This 
is clearly due to lacking metallicity constraints, and because CI 
is considerably more metal-rich than the Sun. If the star had so- 
lar metallicity, the results would be comparable between each 
set. For example, for C3 which has solar metallicity we find a 
comparable radius using the photometric data (1.89 R©) and the 
spectroscopic data (1.85 - 1.93 Rq). For C2 which has sub-solar 
metallicity we find the opposite trend to CI; the fitted radius 
is 2.18 Rq without a metallicity constraint, while adopting the 
lower metallicity values yields smaller radii of 2.11 and 2.13 
Rq. This same trend is found using the results from the other 
pipelines. 

While the absence of a metallicity constraint in general in- 
creases the uncertainties e.g. from 0.05 to 0.11 R© for CI, with 
similar values found for C2 and C3, its absence will also bias the 
final fitted value of radius and hence mass, although still within 
its given uncertainty. For C4 the fitted radius varies by about 2% 
correlating directly with the input T^b but the statistical uncer- 
tainty is 4 - 5%. In this case the absence of a metallicity mea- 
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Table 10. Input spectroscopic observations (analysed along with 
the seismic data) for star CI and the SEEK determination of 
mass and radius. 



Set* 




logs 


[Fe/H] 


R 


0"R 


M 


0"M 




(K) 


(dex) 


(dex) 


(Rg) 




(Mo) 




Si 


5717 


3.96 


+0.35 


2.25 


0.03 


1.45 


0.10 


S2 


5445 


3.84 


+0.13 


2.20 


0.05 


1.33 


0.10 


S3 


5580 


3.81 


+0.19 


2.23 


0.04 


1.37 


0.11 


S4 


5650 


4.10 


+0.36 


2.25 


0.04 


1.44 


0.10 


S5 


5700 


4.10 


+0.13 


2.23 


0.06 


1.39 


0.13 


S6 


5660 






2.13 


0.11 


1.21 


0.20 



surement is unimportant, due to the error in {Av) being a factor 
of 10 larger than for the other stars. 

5.2. Different descriptions of tine input piiysics 

We reported the fitted properties (mass, radius, and age) and 
uncertainties for stars CI - C5 in Table [8] using the reference 
method SEEK and the mean values obtained by combining the 
results from the four or five pipelines. We analysed the data us- 
ing five grid-based seismic methods to test for the reliability of 
the results. Because the five methods are based on different sets 
of physics and input parameters, we can consider the scatter of 
the results from these pipelines as indicative of the systematic 
error that we can expect. Testing all of the possible sources of 
systematic errors in the most coiTect manner would imply gen- 
erating numerous new grids with every combination of physics, 
stellar evolution codes, and parameters, and using the same anal- 
ysis method that SEEK uses. While this is beyond the scope of 
this paper, and beyond the idea behind using grid-based meth- 
ods, we can, however, investigate systematics by varying a few 
of the most important physical ingredients of the models, for 
example, the convective core overshoot parameter, and the in- 
clusion/exclusion of diffusion. We keep in mind that the SEEK 
method covers almost the full range of possible mass, age, ini- 
tial metal fraction, and mixing-length parameter for these stars, 
and the other pipelines use different combinations of EOS, opac- 
ities, nuclear reaction rates, and diffusion of elements (for some 
codes), and these should cover most of the systematics that one 
would expect to find. 

We begun with a stellar model which is described by one set 
of stellar parameters P. SEEK provides central parameters de- 
fined by their distribution in the plane, so we used the SEEK 
parameters as a starting point to determine a single best-fitting 
model using the Levenberg-Marquardt minimisation algorithm, 
and the same input physics as SEEK. Because we have only five 
observations (Teff, iogg, [Fe/H], (Av), Vmax) and in principle five 
fitting parameters (M, t, Zj, Xi or Yi, a), we decided to fix the 
initial hydrogen mass fraction, which in effect allows Yi to vary 
slightly, such that Yj has near solar value, 7, - 0.273 - 0.278 
dSerenelli & BasuL l2010h . Table [TT] lists the values of the best 
fitted parameters using the original input physics in SEEK and 
then for cases la - Id and 2a - 2d, which are described below. 

Once we have found P we change the physical description 
of the model by including convective core overshoot and setting 
its parameter to aov - 0.25, and we search again for a new set 
of parameters that describe the observations best. We do this in 
various ways using the same minimisation algorithm. The first 
time we fit the same four parameters (case la), the second time 
we fit T, Zi and a (case lb), then we fit M, Z\ and a (case Ic), 



and finally we fit the parameters M and a only (case Id). If we fit 
both the mass and age together (case la), then we will determine 
a good model, but it is not possible to test if the new fitted mass 
and age are due to the 'correlation' term between the parameters, 
or due to the new description of the model. For this reason we 
fix the mass and search for a set of parameters including the 
age that adequately fit the observations (case lb) and vice versa 
(case Ic). This yields an estimate of the systematic error on the 
age/mass parameter for a fixed mass/age. We al so include case 

i case I d) to eliminate the M-Y\ correlation e.g. iMetcalfe et al.l 
2009h : lOzel et al.l (1201 Ih . We repeat the same exercise while 
also including He diffusion as described by Michaud & Proffit 
(1993). These cases are denoted by 2a - 2d. 

Inspecting Table fTTl we find a maximum difference of 5% in 
mass and 1 .8% in radius (case 2d), while the largest difference in 
the fitted age is found for case (2b) resulting in a 10% difference 
from the original fitted value. These systematic errors are much 
smaller than the uncertainties given by SEEK in Table |8] While 
we do not claim that these values are typical values for all com- 
binations of physical descriptions in the models, these results 
indicate that the uncertainties are realistic. We may find larger 
differences using the more evolved stars, but the uncertainties 
for these stars are also generally larger. 

6. Comparison of five grid-based approaches 

Using the five Kepler stars and the Sun as test cases for grid- 
based analyses we have shown the following: 

6.7. Surface gravity 

If the seismic data alone are used to determine log g using stel- 
lar models, then we can expect to find systematic differences 
between the pipelines due to the different combinations of in- 
put physics. The difference between all of the results is at most 
0.07 dex, which still provides a very strong constraint on \ogg. 
Combining the results from all of the pipelines yields log g val- 
ues in better agreement with those calculated directly from the 
seismic data (without models) when we make a modest assump- 
tion about Teff. When both atmospheric and seismic constraints 
are available, all of the pipelines yield consistent results for 
\ogg, fitting to within at most 2cr (or 0.05 dex) of either the 
SEEK result, the mean value from all the pipelines, or the model- 
independent value. 

6.2. Radius and mass 

The determination of mass, M, depends mostly on the value of 
the fitted radius, R, because the seismic data constrains the M-R 
correlation. For this reason, similar trends are found for both R 
and M. The systematic error (cr^ys), defined here as the maximum 
difference in fitted values between SEEK and the other pipelines, 
was found to agree within ~2cr. We also found that in general 
for the radius osys increases with the evolutionary state of the 
star, i.e. 2% for the Sun which is mid-MS, and 10% for the most 
evolved star in this study (C4). 

We adopted SEEK as the reference method, which has been 
validated with independent measurements of nearby stars, and 
we find that for each star, the mean parameter value obtained by 
combining the results from all of the pipelines agrees to within 
Icr of the SEEK results, e.g. ^seek - (R) < lo" for all stars. 
This seems to imply that adopting an average value of the stellar 
parameters from several grid-based methods may be the optimal 
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Table 11. The best-fitting parameters for star CI found by using a minimisation algorithm with SEEK physics with some changes 
in the model. 



Descript 


ion 


M 


R 


T 


Y, 


Zi 


a 


xi 






(Mo) 


(Ro) 


(Gyr) 














1.359 


2.253 


3.93 


0.2706 


0.0294 


1.21 


1.26 


(la) ttov 


= 0.25 


1.386 


2.275 


3.81 


0.2745 


0.0255 


1.33 


1.73 


(lb) ttov 


= 0.25 




2.248 


3.70 


0.2774 


0.0226 


1.28 


1.89 


(Ic) Q-ov 


= 0.25 


1.386 


2.264 




0.2739 


0.0261 


1.39 


1.49 


(Id) Q-ov 


= 0.25 


1.417 


2.290 








1.35 


1.18 


(2a) ttov 


= 0.25, He diff 


1.379 


2.265 


3.74 


0.2746 


0.0254 


1.30 


1.75 


(2b) ttov 


= 0.25, He diff 




2.250 


3.55 


0.2763 


0.0237 


1.18 


0.98 


(2c) Q-ov 


= 0.25, He diff 


1.379 


2.271 




0.2745 


0.0255 


1.40 


1.22 


(2d) Q-ov 


= 0.25, He diff 


1.429 


2.295 








1.53 


0.80 



Table 12. Photometric data for CI - C5 from the literature. 



Star 


V 






E(S - VY 




(mag) 


(mag) 


(mag) 


(mag) 


CI 




10.942 


10.901 


0.028 


C2 


10.871 


10.840 


10.959 


0.056 


C3 


11.684 


11.635 


12.004 


0.070 


C4 


11.982 


11.885 


12.160 


0.075 


C5 






12.049 


0.062 



" lUrban et al] fT99^ '' iKharchenkol flOOTh iDroege et alj flOO^ '' KIC 



value to use. Nevertheless, all of the grids provide radius and 
mass values consistent within 2 Act. 

6.3. Age 

In Fig. |4]we showed the fitted mass versus age for stars CI and 
C3 using all of the methods for CI and four methods for C3. 
While the returned value of the age depends upon the description 
of the physics and the mass of the star, in general we found the 
differences among the derived age is due to the mass-age correla- 
tion. Both the RADIUS and Yale-Birmingham pipelines provide 
the smallest uncertainties for CI. These are probably underes- 
timated if one is to consider all possible sources of systematic 
errors. 

When {Sv()2} is available it is optimal to use the methods 
based on oscillation frequencies (SEEK, CESAM2k/Mumbai) 
in order to decrease the uncertainty in the age, especially for 
the less evolved MS stars e.g. for the Sun we find 9.2 + 3.8 Gyr 
without {6vq2) and 4.9 + 0.5 Gyr including it. However, the ac- 
tual fitted value of age does not change very much for the more 
evolved MS stars analysed here. 

7. Additional stellar properties from complementary 
data 

7.1. Determination of distance 

The luminosities corresponding to the models from Table[8]yield 
values of the absolute bolometric magnitude M\,o\. To ob tain the 
absolute V magnitude My we interpolate the tables from lFloweil 
(119961) for the values of T^ff given by VWA in Table [3] (and 
the photometric value for C5) to obtain the bolometric correc- 
tion BC, and apply this coiTection to Mboi. In Table fT2l we list 
the apparent V magnitudes from various sources. To account for 
reddening we use the standard extinction law Ay = 3.1xE(B - 



V), where 3.1 is a typical value dSavaee & MathisL Il979l) and 
E{B - V) are obtained from the KIC. We note that the adopted 
E(jj- V) should be used with caution dMolenda-Zakowicz et al.L 
120091) . The distance modulus is calculated from My and the de- 
reddened V, to give the asteroseismic distance d in parsecs using 
d = iqO-2(v-Mv)+i ^ jjj Table[T3]we list the model L My, the de- 
redden ed V magnitude adopting the values from iDroege et al.l 
(l2006h Vdei, and the distance do^, with the uncertainties arising 
from the luminosity uncertainty. The column with heading c/qi 
shows the distance using V from Kharchenko (2001!). We are 
able to determine distances to these stars with a precision of less 
than 10%. However, for C2 and C4 we find a discrepancy of 14% 
and 12% due to the different reported magnitudes. 

We can also calculat e the distance using o ne of the surface 
brightness relations from lKervella et al] (|2004|) (Table 5), if we 
know the radius. We used the relationship for Teff and V; log = 
3.0415(log reff)2-25.4696 log reft+53.7010-0.2 Vri,r, where 6> is 
the an gular diameter in milliarcseconds. Using the lDroege et al.l 
(l2006h magnitudes, the fitted radii, and the observed T^g, we 
calculated the distance according to this relation c/sb (also given 
in Table [T3ll. and these were found to agree with those using the 
distance modulus. 

7.2. Rotational period and inclination 

The value of v sin / can be determined from the spectroscopic 
analysis (see Table [3]l. Combining this with R allows us to con- 
strain the stellar rotational period Prot - inclination / relation. 
We use the v sin ; values from ROTFIT to determine this rela- 
tionship. With no observational constraints on ; for any star and 
having uncertainties in v sin / of the same order as its value for 
C3 (implying possibly vsini ~ kms"'), the lower bound on 
Prot is unconstrained for all of the stars, while the upper bound 
is poorly restricted for C3. We can place an upper bound on Prot 
for CI, C2, and C4 of 384 64, a nd 1 68 days, respe c tively 

Both lCampante et al.l (l20Ilh and iMathur et al.l (l20Ilh anal- 
ysed the low range of the frequency spectrum to look for signa- 
tures of a rotation period. They estimate Prot for CI, C2, and C3 
of ~ 36, 23, and 27 days, respectively, consistent with our results. 
Adopting these values as Prot constrains the stellar inclination 
angle / = 20°+/5^ for CI, i = 44°^^^ for C2, and / = 25°^^^ for 
C3. 

7.3. Lithium content and age estimate 

We noted that two of the investigated stars, namely CI and C3, 
clearly display an Li i /16707.8 A photospheric absorption line in 
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Table 13. Model luminosities L, absolute V magnitudes My, de-reddened magnitudes Vdei, and asteroseismic distances from the 
standard distance modulus equation (c/o6,oi) and surface brightness relation t/sB- 



Star L(trL) My Vdc, doe ^oi <^sb 
(Lo) (mag) (mag) (pc) (pc) (pc) 

"Cl 4.2(1.1) 3300 10814 SWf^ 324^7^ 318± 20 

C2 5.3 (1.1) 2.967 10.785 366+| 347!^ 371± 24 



C3 3.6(1.2) 3.427 11.787 470+^^ 396™ 475± 31 
C4 20.0(1.1) 1.510 11.927 1212^^3 1068+]^ 1231±95 
C5 7.8 (1.1) 2.526 11.857 735+g 741± 61 



the FIES spectra (Fig. |5]l, while for the remaining two stars the 
lithium line is not detectable. Lithium is burned at relatively low 
temperatures in stellar interiors (~ 2.5 x 10^ K). As a conse- 
quence, it is progressively depleted from the stellar atmospheres 
of late-type stars when mixing mechanisms pull it deep into their 
convective layers. Therefore, its abu ndance can be used for es- 
timating the stellar age, as shown by ISestito & Randic hi (120051) 
for stars belonging to 22 open clusters spanning the age range 
0.005-8 Gyr. 

We measured the equivalent width of the lithium line, cor- 
recting for the small contribution of the nea rby Fe-I/16707.4A 
line as suggested by ISoderblom et al.l (Il993h . and found Wu - 
102 + 10mA and Wu = 55 + 10mA for Cl and C3, respec- 
tively. For C2 and C4, we estimated an upper limit for the lithium 
equivalent width Wu < 10 mA as the product of the error of 
the normalised flux per spectral point and the integration width 
(- 1 A). 

We derived a lithium abundance logn(Li) - 2.6 + 0.1 
and logn(Li) = 2.4 + 0.1 for Cl and C3 by interpolation of 
flie NL TE curves of growth tabulated by iPavlenko & Magazziil 
(119961) . where by definition logn(H) = 12. For C2 and C4, the 
lithium abundance is log n(Li) < 1 .9. The correlation of lith ium 
abundance and age established bv ISestito & RandichI (120051 see 
their Fig. 7), suggests an age of 0.1 - 0.4 Gyr and 1-3 Gyr for 
Cl and C3, respectively, while their Table 3 suggests an age for 
C2 and C4 corresponding to the most evolved clusters in their 
study (M67 at 5 Gyr). The asteroseismic ages of 3 Gyr for C2 
and 1 Gyr for C4 imply at least evolved MS stars for masses of 
1.4 and 1.9 M©, respectively. 

New Wu measurements for solar-like stars in galactic open 
clusters in the age range 1-8 Gyr, have shown that the low solar 
lithiu m abundance is not the stand ard for a star of that a ge and 
mass (lRandichil2010h . In particular. iPasquini et al.l (l2008l) found 
a spread of lithium abundance in stars of the solar-age cluster 
M 67, ranging from ~0.4 to ~2.0 dex. Randich (2010) also found 
that the average abundance for various clusters of similar age 
varies. For stars in the Tetj range 5750-6050 K, some evidence 
of bimodality in lithium depletion after 1 Gyr seems to emerge 
from these new data, with some clusters following the solar be- 
haviour of abundance decay and other ones forming a plateau 
with log«(Li) = 2.2 - 2.4. It is now unclear what the driving 
mechanism for depletion is, but it appears that some other pa- 
rameter apart from mass and age is playing a role. 

This could reconcile the lithium abundance log«(Li) = 
2.4 + 0.1 of C3 with the higher age derived by asteroseismology. 
However, for Cl, whose temperature is significantly lower (more 
lithium depletion), the lithium abundance of log n(Li) = 2.6±0.1 
is not compatible with the age of 4 Gyr deduced by asteroseis- 
mology. We note that Cl is the most metal -rich star studied in 
this paper. 
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Fig. 5. A region of the optical spectra which contains the Lithium 
/16707.8A line. From top to bottom we show the spectra from 
Cl, C3 and a high-resolution solar spectrum (Ganymede taken 
in 2007 with HARPS). 



8. Summary and Conclusions 

In this work we analysed the five solar-type stars KIC 1 1395018, 
KIC 10273246, KIC 10920273, KIC 10339342, and 
KIC 11234888, referred to as Cl - C5, respectively, based 
on more than eight months of short-cadence Kepler data. The 
global seismic quantities ((Av) and v^ax) coupled with atmo- 
spheric parameters derived from photometry and spectra taken 
with the NOT telescope yield stellar properties with an average 
precision of 2% in mean density, ~0.03 dex in surface gravity, 
2-5% in radius, 7-11% in mass, and 35% in age (Table [T4li. 

We used five grid-based approaches based on stellar mod- 
els to estimate systematic errors, since the grids are all based 
on different physics and evolution codes. We found an agree- 
ment between all of the methods to within 2.4cr for mass and 
radius, and these values were often comparable to the statisti- 
cal uncertainties (see Tables |5] and |6]l. However, we also tested 
some specific sources of systematic errors arising from includ- 
ing convective core overshoot and settling of helium. We found 
that the fitted stellar parameters (radius, mass, and age) changed 
by an amount smaller than the statistical uncertainties reported 
in Table[T4l These changes were also smaller than the dispersion 
among the results from the five grid-based methods, and we can 
therefore conclude that using several different methods provides 
an estimate of the systematic errors. However, since the value 
of the age is highly model-dependent, and we only tested two 
physical ingredients in the models, the systematic errors in the 
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age could be larger. In Table [TSl we summarise the estimates of 
the systematic errors osys from various sources. 

The spectroscopic data were analysed by five groups inde- 
pendently, each providing atmospheric parameters that varied by 
more than 200 K in T^s and 0.2 dex for [Fe/H]. The differences 
obtained can be attributed to the quality of the spectroscopic data 
(low S/N and intermediate resolution). However, we found that 
the results are in general correlated (Fig.O, indicating that with 
higher resolution and higher S/N data, the results should be in 
much better agreement. The different results obtained, however, 
also allowed us to study the influence of these atmospheric pa- 
rameters for determining the radius and mass of the stars. We 
found that the fitted radius depends on both the adopted T^ff and 
[Fe/H], and the different values led to a discrepancy of about 
2.5% in stellar radius and 8% in stellar mass. In the absence 
of an [Fe/H] measurement we found biased results in the fitted 
radius and mass for stars with non-solar metallicity (metal-rich 
CI and metal-poor C2), and larger uncertainties. Higher-quality 
spectra should be obtained for these targets in the near future to 
match the exceedingly high quality of the Kepler data. 

We investigated the role of including (5v()2) as an observa- 
tional constraint with (Av) and Vmax (Sect. |43] and Table |9l), and 
we found that (5vo2) is important primarily for the uncertainty in 
age for CI, C5, and the Sun. For C5 this is due to the absence of 
an [Fe/H] measurement. We suspect that its larger impact for the 
Sun, is because it is in the middle of its MS lifetime where the 
other atmospheric observables change relatively slowly with the 
age, thus providing weaker constraints. 

Coupling photometric magnitudes with the model luminosi- 
ties, we derived distances with a precision of less than 10%, and 
for C3 and C4, systematic errors of 14% and 12% respectively 
arise from the V magnitudes reported by different authors. Later 
on in the mission, parallaxes may be obtained for the unsaturated 
stars. These independent determinations of distances should help 
to resolve the differences in the V magnitudes and/or reddening. 

By coupling the derived radius with the observed v sin /, lim- 
its can be placed on the rotational period f rot of the stars. For 
CI, C2, and C4 we can impose an upper bound on Prot of 386, 
67, and 177 days, respectively. Rotational period, inclination, 
and rotational velocity can all be determined in an independent 
manner by studying either the stellar spot distribution using the 
time series data, the frequency splittings in the power spectrum, 
or by stud ying the low frequency range of the p ower spectrum. 
ICampante et al. (.2011) and Mathur et al. (i201lh study the low 
frequency range of the power spectra and estimate rotational pe- 
riods of 36, 23, and 27 days for CI, C2, and C3, respectively. 
Adopting these as Prot constrains the stellar inclination angle 
i = 20°+/^^ for CI, / = 44°+4^ for C2, and ; = 25°+,^,^ for C3. 

The amount of lithium absorption in the atmospheric spectra 
can be used to estimate the age of the star. For the stars C2 and 
C4 no Li absorption is seen, indicating evolved MS stars, and 
this is confirmed by the seismic analysis. For C3, the Li abun- 
dance indicates an age that is inconsistent with the seismic value, 
but these could be reconciled, by consideri ng new results w hich 
show a bimodal distribution of Li depletion (lRandichl.l2010l) . For 
the metal-rich star CI, however, the relatively high Li abundance 
indicates a non-evolved MS star, while the asteroseismic analy- 
sis suggests an incompatib le v alue of ~ 4 Gy r. 

ICampante et alJ(l201 ih and lMathur et al.1(l201 ih analysed the 
time series and presented the oscillation frequencies for CI, C2, 
C3, and C5. Using the individual frequencies will yield more 
precise stellar properties, as well as studies of the stellar in- 
terior (Brandao et al . 2012, Dogan et al. 2012) as shown by 
iMetcalfe et al.l (l2010h . Up to now, such precise values have only 



been possible for members of detached eclipsing binaries, the 
Sun, and a few very bright stars. Thanks to the significant im- 
provement in data quality from Kepler, this is about to change. 
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Appendix A: Oscillation frequencies 

In Table lA. II we list th e p ublished oscil l ation frequencies from 
ICampante e'taD(l20 lib and lMathur et al.l (1201 ih for stars CI, C2, 
C3, and C5. 

Appendix B: Atmospheric analysis methods 

B.1. BIA 

The MOOG code (version 2002: ISnederilll973h determines the 
iron abundance under the assumption of local thermodynamic 
equilibrium (L TE), using a grid of ID model atmospheres by 
lKuruczl([T993h . The LTE iron abundance of C 1 - C4 was derived 
from the equivalent widths of 55-68 Fe-I and 1 1 Fe-II lines in the 
4830-68 10 A range, measured with a Gaussian fitting procedure 
in the I RAlt^task splot. For the analysis, we used the same list of 
lines as lBiazzo et al.l (1201 ll) and we followed their prescriptions. 
In particular, Teff and were determined by requiring that the 
iron abundance be independent of the excitation potentials and 
the equivalent widths of Fe-I lines. The surface gravity was de- 
termined by requiring ionization equilibrium between Fe-I and 
Fe-II. The initial values for Teff, log g and were chosen to be 
solar (reif = 5770 K, logg = 4.44 dex, and|t = 1-1 km s^'). 

B.2. ROTFIT 

The IDlB code ROTFIT dFrasca et all l2006l) performs a simul- 
taneous determination of Tefr, log g, and [Fe/H] for a star, as well 

^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of the Universities 
for Research in Astronomy, inc. (AURA) under cooperative agreement 
with the National Science Foundation. 

* IDL (Interactive Data Language) is a registered trademark of ITT 
Visual Information Solutions. 
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Table 14. Summary of stellar properties for KIC 11395018 (CI), KIC 10273246 (C2), KIC 10920273 (C3), KIC 10339342 (C4), 
and KIC 1 1234888 (C5), derived from atmospheric and mean seismic parameters, and stellar models. 





CI 


C2 


C3 


C4 


C5 


»,(kgm-') 


175±2 


185±1 


254±1 


39±16 


135±1 


''<p>MR (kg m-3) 


174±4 


189±2 


257±5 


46±4 


140±2 


"loggv (dex) 


3.86±0.03 


3.88±0.03 


3.94±0.03 


3.47±0.03 


3.79±0.04 


''log gMR (dex) 


3.88+0.02 


3.88±0.02 


3.97±0.03 


3.52±0.04 


3.82±0.03 




2.23±0.04 


2.11±0.05 


1.90±0.05 


3.81±0.19 


2.44±0.14 


M (Mo) 


1.3 /±0. 1 i 


i.Z0±U. to 


i.ZD±U.i3 


1 net 1 n 1 o 
1. /9±U. 12 


i.44±U.Z0 


T (Gyr) 


3.9±1.4 


3.7±0.7 


4.5±1.8 


1.1±0.2 


2.6+0.9 


^T(ey} (Gyr) 


4.5±0.5 


3.7±0.6 


5.0±1.9 




1.6+0.2 


L(Lo) 


4.2±1.1 


5.3±1.1 


3.6±1.2 


20.0±1.1 


7.8+1.1 


reffmodcl (K) 


5547 


6047 


5789 


6255 


6180 


in 




4.4+46 


25+65 






-PROTmax (days) 


384 


64 




168 




''/'ROTesi (days) 


36 (6) 


23 


27 




19-27 (5) 


rf(pc) 


318+3^ 


366+'* 


470!- 


1212^33 





Subscripts v and MR indicate that the value was obtained directly from the data and from the models, respectively, 
''r^av) is when {Svq2) is included as an observa tiona l constraint. In this c ase the uncertainties in logg, R, and M for C5 reduce by a factor of three. 
''PaoTsa as reported bvlCamp ante et al] ( l201lh and lMathur et al.l i201 ih . with uncertainties in parenthesis. 

Table 15. Estimates of systematic errors in the stellar properties, given in CGS units for g, solar units for R and M, and Gyr for age, 
with % values shown in parentheses. 





CI 


C2 


C3 


C4 


C5 


""logs.grid 


0.02 


0.05 


0.01 


0.04 


0.02 


C" R.grid 


0.06 (3) 


0.22(11) 


0.07 (4) 


0.38(10) 


0.12(5) 


^^M.grid 


0.09 (7) 


0.43 (34) 


0.06 (5) 


0.09 (5) 


0.18(12) 


^r.grid 


1.4 (36) 


1.7 (38) 


1.8 (40) 


0.2(18) 


0.4(16) 


0" R.spec 


0.05 (2) 


0.08 (4) 


0.03 (2) 


0.18(5) 






0.13 (9) 


0.16(13) 


0.03 (2) 


0.12(7) 




O" R.phys 


0.06 (3) 










C M.phys 


0.05 (4) 










^r.phys 


0.35 (9) 











The subscripts 'grid', 'spec', and 'phys' mean estimates of systematic errors from using different grids/pipelines, different atmospheric constraints, 
and different descriptions of the physics in the models, respectively (see Sects. |4]and|5]for details). 



as its projected rotational velocity v sin /, by comparing the ob- 
serv ed spectrum with a library of spectra of reference stars (see, 
e.g.. lKatz et 10119981: ISoubiran et anil998h . We estimated the 
stellar parameters and vsin; as the mean values for the 10 ref- 
erence stars whose spectra, artificially broadened with v sin / in 
steps of 0.5 kms"', most closely resembled the target spectrum, 
and each is quantified by a x' measure. We adopted their stan- 
dard deviation as a measure of the uncertainty. We applied the 
ROTFIT code to all echelle orders which cover the range 4300- 
6680 A in the observed spectrum. The adopted estimates for the 
stellar parameters come from a weighted mean of the values for 
all the individual orders, where more weight is assigned to the 
best fitted or higher S/N orders, and the "amount of information" 
contained in each spectral region expressed by the total spectral- 
line absorption. The standard error of the weighted mean was 
adopted as the uncertainty estimate for the final values of the 
stellar parameters. With the same code, we also performed an 
MK spectral classification for our targets. 

B.3. NIEM 

The NIEM analysis follows the methodology presented in 
iNiemczura. Morel & AertsI (l2009h and relies on an eflicient 
spectral synthesis based on a least squares optimisation algo- 



rithm (see lTake"dalll995h . This method allows for the simulta- 
neous determination of various parameters involved with stellar 
spectra and consists of the minimisation of the deviation between 
the theoretical flux distribution and the observed normalised one. 
The synthetic spectrum depends on the stellar parameters: ef- 
fective temperature Tefj, surface gravity logg, microturbulence 
^t, rotational velocity v sin /, radial velocity vr, and the relative 
abundances of the elements. The first three parameters are ob- 
tained before the determination of abundances of chemical el- 
ements and they are considered as input parameters. All the 
other above-mentioned parameters can be determined simulta- 
neously because they produce detectable and different spectral 
signatures. 

Tetf , log g, and are adjusted by comparing the abundances 
determined from the unblended Fe-I and Fe-II lines. We require 
that the abundances measured from Fe-I and Fe-II yield the same 
result. The absorption lines of neutral iron depend only on Tefj, 
^t, and iron abundance, and are practically independent of the 
surface gravity. On the other hand the lines of ionised iron are 
sensitive to the temperature, metallicity and most of all to grav- 
ity. First, we adjust until there is no coiTelation between iron 
abundances and line intensity for the Fe-I lines. Second, Tefj is 
adjusted until there is no trend in the abundance versus excita- 
tion potential of the atomic level causing the Fe-I lines. Both 
and reft are not independent, however the dependence on is 
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Table A.l. PubUshed frequencies given in ^Hz for stars KI C 11395018 (CI), K IC 10273246 (C2), KIC 10920273 (C3), and 
KIC 11234888 (C5), as reported bv lcampante et all (1201 ih and lMathur et"aLl(l201lh . 



/ 


KIC 11395018 


KIC 10273246 


KIC 10920273 


KIC 11234888 




CI 




C2 




C3 




C5 









686.66 ± 


0.32 


737.90 ± 


0.30 


826.66° H 


tO.25 


627 


67 + 


0.19 





732.37 ± 


0.18 


785.40 ± 


0.20 


882.77 + 


0.20 


669 


35 + 


0.16 





779.54 ± 


0.14 


833.90 ± 


0.20 


939.58 + 


0.16 


711 


63 + 


0.15 





827.55 ± 


0.15 


883.50 ± 


0.20 


997.14 + 


0.18 


753 


64 + 


0.20 





875.40 ± 


0.16 


932.70 ± 


0.50 


1054.33 : 


±0.30 


794 


56 + 


0.20 





923.16 ± 


0.19 


981.10 ± 


0.30 


1111.51 : 


± 0.25 


836 


83 + 


0.21 





971.05 ± 


0.28 


1030.70 


± 0.40 


1170.77° 


+ 0.33 


877 


80° : 


± 0.22 









1079.30 


± 0.20 


1226.34° 


+ 0.33 








1 


667.05° : 


t 0.22 


622.80 ± 


0.20 


794.65* H 


t 0.32 


506 


72 + 


0.20 


1 


707.66 ± 


0.19 


661.90 + 


0.50 


838.61* H 


tO.25 


563 


30 + 


0.14 


1 


740.29* : 


t 0.17 


695.75* : 


t 0.27 


914.52 + 


0.16 


594 


83 + 


0.16 




763.99 ± 


0.18 


724.70 + 


0.20 


968.19 + 


0.13 


686 


34 + 


0.17 


1 


805.74 ± 


0.13 


764.30 + 


0.30 


1023.58 : 


± 0.14 


741 


10 + 


0.18 


1 


851.37 + 


0.11 


809.80 + 


0.20 


1079.10: 


± 0.31 


815 


43 + 


0.21 


1 


897.50 ± 


0.15 


857.30 + 


0.20 


1135.36" 


+ 0.31 


855 


67 + 


0.21 


1 


940.50 ± 


0.15 


905.60 + 


0.30 












1 


997.91 ± 


0.33 


950.00 + 


0.30 












1 






1008.60 


tO.40 












1 






1056.30 


±0.20 












1 






1103.30 


±0.40 












2 


631.19° : 


t 1.36 


688.50 + 


0.70 


822.39° H 


t 0.28 


582 


84° : 


± 0.21 


2 


680.88 ± 


0.45 


734.80 + 


0.60 


873.10°'' 


+ 0.32 


624 


65 + 


0.18 


2 


727.78 ± 


0.30 


779.50 + 


0.40 


934.49 + 


0.22 


708 


67 + 


0.19 


2 


774.92± 


0.16 


830.30 + 


0.40 


992.44 + 


0.13 


751 


84 + 


0.21 


2 


823.50 ± 


0.16 


880.60 + 


0.50 


1049.36 : 


±0.39 








2 


871.29 ± 


0.21 


927.50 + 


0.40 


1106.76 : 


±0.34 








2 


918.10 + 


0.28 


977.60 + 


0.40 












2 


965.83 ± 


0.23 


1025.30 


± 1.30 












2 


1016.61° 


f ± 0.73 


1073.70 


± 0.20 













2 ... 1122.70°'' + 0.40 



The frequencies were fitted with 11, 10, 10, and 1 1 different methods, and a criteria was imposed to obtain a minimal frequency set representative 
of the results from most of the fitters. 
° Frequencies not among the minimal set. 
''1=1 mixed mode. 

" Mode close to the second harmonic of the inverse of the long-cadence period. 
'' Possible / = 2 mixed mode introduced a posteriori. 
Uncertain identification. 

Mode tagging changed a posteriori from / = to / = 2. 



stronger, so we adjust this parameter first. The gravity is obtained 
by fitting Fe-II and Fe-I lines and requiring the same abundances 
from both neutral and ionised lines. 

B.4. SOU 

The SOU parameters were derived starting with the automatic 
measur ement of equivalent widths of Fe-I and Fe-II lines with 
ARES (ISousa et al 1 I2007h and then imposing excitation and 
ionization equilibrium using a spectroscopic analysis in LTE 
with the help of the code MOOG (Sneden, 1973) and a GRID 
of Ku rucz Atlas 9 plane-parallel model atmospheres dKuruczl 
Il993h . The Fe I and Fe II line list is composed of more than 
300 lines that were individually tested in high resolution spectra 
to check its stability to an automatic measurement with ARES. 
The atomic data of th e lines were obtain ed from the Vienna 
Atomic Line Database (iKupka et al l ll999h butthe log gf were 
recomputed through an inverse analysis of the solar spectrum 



allowing in this way to perform a differential analysis relative 
to the Su n. A full description of the method can be found in 
ISousa et a l. (2008). The errors on the SOU parameters are ob- 
tained by quadratically adding 60 K, 0.1 and 0.04 dex to the 
internal errors on Tss, iogg, and [Fe/H], respectively. These val- 
ues were obtained considering the typical dispersion plotted i n 
each comparison of parameters presented in lSousa et alj (l2008h . 
A more complete discussion about the errors derived for this 
spectroscopic method can be seen in lSousa et al. l (l20Tll) . 

B.5. VWA 

The VWA software dBruntt et al.L 1201 Ol) relies on spectral syn- 
thesis of the region around individual Fe lines, but includes the 
contribution from weakly blending lines. The selection of lines 
was chosen from the degree of blending. These were then fitted 
with a sy nthetic spectrum. We ad opted the MARCS model atmo- 
spheres (iGustafsson et al.L l2008l) and used the atomic line data 



16 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 



from the Vienna Atomic Line Database dKupka et a ll Il999h . 
Each line fit was inspected in great detail and bad fits discarded, 
resulting in between 100-150 Fe-I/Fe-II lines that were used 
in the parameter determination. The oscilla tors strengths log gf 
were corrected relative to the Sun as done in lBruntt et al.l (l2010l) . 
The values quoted in the Kepler Input Catalogue were used as 
initial guesses for the parameters of the atmosphere. The param- 
eters were then refined through several iterations where each pa- 
rameter was adjusted so that any correlations between the abun- 
dances of Fe-I and equivalent width (EW) and excitation poten- 
tial (EP) were removed. Also, we required agreement between 
the Fe-I and Fe-II abundances. Only lines with an EW < 90 mA 
were used, since these lines are the most sensitive to changes in 
the fundamental parameters of the star. 

To calculate the uncertainties of Teff , log g, and the model 
parameters were changed one at a time, until at least a 3cr devia- 
tion was produced on the slope of the Fe-I-abundance vs. EP or 
EW, or between the Fe-I/Fe-II abundances. From these, the Icr 
error was calculated, giv ing the internal prec ision of the parame- 
ters calculated in VWA dBruntt et al.Ll2008l) . The uncertainty on 
[Fe/H] was calculated as the standard error on the mean from Fe- 
I lines and adding quadratically the uncertainty from the model 
Teff, \ogg, and^f 

To put further constraints on the \ogg value of the targets 
stars, we also used individual fits of the pressure-sensitive Mg ib 
lines and the Ca Unes at .i6122A, .i6162A and ^16439 A. Each 
line was fitted with three different synthetic spectra, each with a 
different value for log g. For each fit the was calculated and 
from this the best value of log g was dete rmined. The metho d 
is described in greater detail in a paper bv iBruntt et al.l (l2010l) . 
The final values for log g were taken as the weighted mean of the 
values found from matching the Fe-I/Fe-II abundances and from 
fitting the pressure-sensitive lines. 
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